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Abstract 

We present a calculation of the effective cross section Ue//, an important ingredient 
in the description of double parton scattering in proton-proton collisions. Our theoretical 
approach makes use of a Light-Front quark model as framework to calculate the double parton 
distribution functions at low-resolution scale. QCD evolution is implemented to reach the 
experimental scale. The obtained CTe//, when averaged over the longitudinal momentum 
fractions of the interacting partons, Xi, is consistent with the present experimental scenario. 
However the result of the complete calculation shows a dependence of Oeff on li, a feature 
not easily seen in the available data, probably because of their low accuracy. Measurements of 
de// in restricted Xi regions are addressed to obtain indications on double parton correlations, 
a novel and interesting aspect of the three dimensional structure of the nucleon. 
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1 Introduction 


Multi Parton Interactions (MPI), occurring when more than one parton scattering takes place 
in the same hadron-hadron collision, have been discussed in the literature since long time ago 
[T] and are presently attracting considerable attention, thanks to the possibilities offered by the 
Large Hadron Collider (LHC) (see Refs. [21 [3l 01 O [6] for recent reports). In particular, the 
cross section for double parton scattering (DPS), the simplest MPI process, depends on pecu¬ 
liar non-perturbative quantities, the double parton distribution functions (dPDFs), describing 
the number density of two partons with given longitudinal momentum fractions and located 
at a given transverse separation in coordinate space. dPDFs are naturally related to parton 
correlations and to the three-dimensional (3D) nucleon structure, as discussed also in the past 

m- 

No data are available for dPDFs and their calculation using non perturbative methods is 
cumbersome. A few model calculations have been performed, to grasp the most relevant features 
of dPDFs [SllHllTO]. In particular, in Ref. m a Light-Front (LF) Poincare covariant approach, 
able to reproduce the essential sum rules of dPDFs, has been described. Although it has not yet 
been possible to extract dPDFs from data, a signature of DPS has been observed and measured 
in several experiments mi Ill da o la 06]: the so called “effective cross section”, o-g//. 
Despite of large errorbars, the present experimental scenario is consistent with the idea that 
(Teff is constant w.r.t. the center-of-mass energy of the collision. 

In this letter we present a predictive study of cTeff which makes use of the LF quark model 
approach to dPDFs developed in Ref. m- 

The definition of Ue// is reviewed in the next section where the present experimental situation 
is also summarized. Then the results of our approach are presented critically discussing the 
dynamical dependence of aeff in view of future experiments. Conclusions are drawn in the last 
section. 


2 The effective cross section 

The effective cross section, erg//, is defined through the so called “pocket formula”, which reads, 
if final states A and B are produced in a DPS process (see, e.g., 0): 


^eff 


,pp pp 


m a A a 


B 


a 


pp 

double 


( 1 ) 


m is a process-dependent combinatorial factor: m = 1 if A and B are identical and m = 2 if they 
are different. differential cross section for the inclusive process pp' —)• A{B) + X, 

naturally defined as: 
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( 2 ) 


{xi,x[,fii) = '^F[{xi,ni)F^ {xr,^ii)afk{xi,xi',^ii) , 
i,k 

(x2,X2,/i2) = '^F^{x2,^i2)Ff {x2',^l2)^fl{x2,X2',^^2) , (3) 

where is a one-body parton distribution function (PDF) with i,j,k,l = {q,q,g}, fJ-i( 2 ) is 
the renormalization scale for the process A{B), the remaining ingredient in Eq. ([T|), 

appears in the natural definition of the cross section for double parton scattering: 




C^Zublei^i^^l’^2,X2,fJ.l,g2) dxidx'^dx2(1x2 , 


(4) 


and reads: 
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.pp 

double 


{xi,x[,X2,X2,Hl,H2) = ^ X] / ^b(^l>^2;k_L,/ri,/r2)<5-)t(xi,Xi)(Tj7(x2,X2) 


X Dki{x[,X2;-'k±,g,i,n2) 


dk i 


(2vr: 


i2 • 


(5) 


In the above equation, k_i_ (—k_i_) is the transverse momentum unbalance of the parton 1 (2), 
conjugated to the relative distance r_i_ (the reader should not confuse k_i_ with the intrinsic mo¬ 
mentum of the parton, argument of transverse momentum dependent parton distributions). The 
quantity Dij{xi,X 2 ', k_i_), called sometimes 2 GPDS [IS [18], is therefore the Fourier transform of 
the so called double distribution function, Zljj(xi,X 2 ;r_i_), which represents the number density 
of partons pairs i,j with longitudinal momentum fractions xi,X 2 , respectively, at a transverse 
separation in coordinate space. dPDFs, describing soft Physics, are nonperturbative quanti¬ 
ties. 

Two main assumptions are usually made for the evaluation of dPDFs: 
a) factorization of the transverse separation and the momentum fraction dependence: 


Dij{xi,X2,'k^,g) = Dij{xi,X2,g)T (k_L, ^); (6) 

b) factorized form also for the xi,X 2 dependence: 


Dij{xi,X2,g.) = Fi{xi,g) Fj{x2,g) 6{l - xi - X2)(l - xi - X 2 T ■ (7) 

The expression 0(1 — xi — X 2 )(l — xi — X 2 )"', where n > 0 is a parameter to be hxed phenomeno¬ 
logically, introduces the natural kinematical constraint xi + X 2 < I (in Eqs. ([6|) and ([7|) the 
same scale g, = g,i, g ,2 is assumed, for brevity). 
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Figure 1; Centre-of-mass energy dependence of fTg// measured by different experiments using 
different processes [HKiailsKllKlSKIS]. The figure is taken from m 


One comment about the physical meaning of aeff is in order. In Eq. ([T]), if the occurrence of 
the process B were not biased somehow by that of the process A, instead of the ratio OB/oeff 
one would read (JBlcTineh representing the probability to have the process B once A has taken 
place assuming rare hard multiple collisions. The difference between (Te// and measures 
therefore correlations between the interacting partons in the colliding proton. 

Let us discuss now the dynamical dependence of Ue// on the fractional momenta xi, x'l, X 2 , X 2 - 
By inserting Eqs. (2-5) in Eq. (1), and omitting the dependence on the renormalization scales 
for simplicity, one gets the following expression for (Jeff- 


^2) 3^2) — 


{(a^i') < 7,4 (3:1 , XU )} { Ylj,i Pj {x2)Ff' {X2' ) afi {x2 , )} 

<74(3^2, x'2) f Ai(xi, X2; k_L) L>fc/(x'i, x'2; -k_L)^ 


( 8 ) 


Eq.dH) clearly shows the dynamical origin of the dependence of on the fractional mo¬ 
menta xi,x(,X 2 ,X 2 . Even within the “zero rapidity region”, (y = 0), where xi = x[, X 2 = X 2 , 
such a dependence, although simplified, is still effective. 

Assuming that heavy flavors are not relevant in the process, the dependence on the “parton 
type”, i = q,q,g, of the elementary cross section is basically [T9] : 


aij{x,x') = Cija{x,x') , 


(9) 
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where <y{x, x') is a universal function, and Cij are color factors which stay in the ratio: 


Cgg : C,g :Cgg = l: (4/9) : (4/9)2 . 


Using Eq. ([9]), Eq. ([8|) simplifies considerably: 


( 10 ) 


*^e// (^1) ) ^2) ^ 2 ) 


Y.i,k,j,l^iixi)Pk{Xi)Fj{x2)Fi{x'2)CikCjl 

CikCjl f Dij{xi, X2; li±)Dki{x^, X2] ~k_i_) 


( 11 ) 


The present experimental scenario is illustrated in Fig.l. The experiments [IIKIIKisKIllia 
IS], at different values of the center-of-mass energy,-^, and with different final states, explore 
different regions of Xj. Experiments at high ^/s access low Xi regions, in general. The old AFS 
data [H] are in the valence region (0.2 < x, < 0.3), the Tevatron data [131114j are in the range 
0.01 < Xi < 0.4 while the recent LHC data [El cover a lower average Xj range and are 
dominated by the glue distribution. 

Remarkably the experimental evidences are compatible with a constant value of (Jeff in 
Eq-dH)) the Xj-dependence being probably hidden within the experimental uncertainties. In fact 
one should stress that the knowledge of the Xj-dependence of cJe// would open the access to 
information on the x^-dependence of the dPDFs T>ij(xi, X 2 ; r_i_), entering the definition of Cg//: 
a direct way to access the 3D nucleon structure [7]. Nowadays, the aspects of the 3D nucleon 
structure related to the transverse position of partons are investigated through hard-exclusive 
electromagnetic processes, such as deeply virtual Compton scattering (DVCS), extracting the 
Generalized Parton Distributions (GPDs) (see Ref. [20] for recent results). The information 
encoded in DPS, dPDFs and in cr^ff, in its full Xj dependence, are anyway different and comple¬ 
mentary to those provided by GPDs in impact parameter space. While the latter quantities are 
one-body densities, depending on the distance, of the interacting parton with given x, from the 
transverse center of the target, in DPS one is sensitive to the relative distance between two par- 
tons with given longitudinal momentum fractions. In other words, the investigation of dPDFs 
from DPS, is relevant to know, for example, the average transverse distance of two fast partons 
or two slow partons: a very interesting dynamical feature, not accessible through GPDs. 


3 Light-Front quark model calculation of the effective cross sec¬ 
tion 

dPDFs have a non-perturbative nature, and, at present, cannot be calculated in QCD. However 
they can be explicitly calculated, at a low resolution scale, Qq ~ Aqcd, using quark models, as 
extensively done for the usual PDFs. The results of these calculations should be then evolved 
using perturbative QGD (pQGD) in order to match data taken at a momentum scale Q > Qq. 
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The procedure is nowadays well established (see, e.g., Ref. [21] and references therein). The 
QCD evolution procedure of dPDFs (from Qq to Q > Qo) is well known |22l[23|, and currently 
implemented in a systematic way (see Ref. laisii and references therein). 

The first model calculations of dPDFs in the valence region, at the hadronic scale Qq, have 
been presented in a bag model framework |S|, and in a constituent quark model (CQM), [Sj. 
Of course CQM have the peculiar advantage of including correlations in a way consistent with 
the quark dynamics, from the very beginning, a property that the bag model cannot fulfill. 

In particular the fully Poincare covariant Light-Front model approach we developed in Ref. 
m respects relevant symmetries, broken in the descriptions of Refs. PE], allowing for a correct 
evaluation of the Mellin moments of the distributions and, consequently, for a precise pQCD 
evolution to high momentum transfer. In this way our model calculations can be relevant for 
the analysis of high-energy data. 

The model, extensively applied to the evaluation of different parton distributions, (see, e.g.. 
Refs. |25l[26l[27] and references therein), is a good candidate to grasp the most relevant features 
of dPDFs. (see Ref. [25] for details). For the present study it is enough to recall that the 
proton state is given by a spatial wave function and an SU(6) symmetric spin-isospin part. 
The spatial part is numerical solution of a relativistic Mass equation, dynamically responsible 
for the presence of correlations between the two quarks in the CQM wave function (a non- 
relativistic version of the model was introduced in Ref. [28]). The Light-Front calculations of 
Djj(xi,X 2 ,k_L,/i), in Ref. [lO], shows that the factorization of Eq. Q is basically valid, but 
the common assumption of Eq. ([7|) is strongly violated. Besides, the strong correlation effects 
present at the scale of the model are still sizable, in the valence region, at the experimental 
scale, i.e. after QCD evolution. At the low values of x, presently studied at the LHC, the 
correlations become less relevant, although their effects are still important for the spin-dependent 
contributions to unpolarized proton scattering. 

We have explicitly calculated single and double parton distributions entering Eq. ffTTTi . and 
then (Teff relying on the natural assumption Eq. (lIOp only. We adhere, in addition, to the 
common choice of a single renormalization scale /ii = = /^O; where has to be interpreted, 

in the present approach, as the hadronic scale where only valence quarks u and d are present. 
Considering the symmetries of our model, one has u{x,fio) = 2d{x,fio), Du,uixi,X 2 ,k±, fio) = 
‘^Du,dixi,X 2 ,k±, Ho) and Eq. (fTT]l simplifies to 


^effixi,Xi,X2,X2,Ho) — 


81u{xi,Ho) u{x[,Ho) u{x2.Ho) u{x'2,Ho) 
64 f Duu(xi, X 2 ;k±, Ho) Duu(x[, x' 2 ; -k±, Ho) 


( 12 ) 


Since the experimental data in the valence region are mainly restricted at zero rapidity 
(y = 0), where = x', one remains with 


^effix\^X 2 ^ Ho) — 


8Iu(xi,/ro)^u(x2./io)^ 
64 / D„„(xi, X 2 ; k_L, 


(13) 
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Figure 2: aeff{xi,X 2 ,Q‘^) for the values of xi,X 2 measured in Ref. [H]. Left panel; hadronic 
scale; right panel: = 250 GeV^. 


In order to illustrate our results we will concentrate on the valence region where the present 
model is more predictive. In particular we concentrate on the kinematics of the old AFS data 
m, which means y = 0 (xi = x'i,X 2 = X 2 ) and 0.2 < xi ,2 < 0.3. The average momentum scale, 
again assumed to be the same for the processes initiated by the two different collisions, turns 
out to be ~ 250 GeV^. The results of the calculations are shown in Fig. 2, at the scale of 
the model, /Tq ~ 0,1 GeV^, and after non-singlet evolution to (details on the hxing of the 
hadronic scale and on the calculation of the QGD evolution can be found in Ref. |10]1. 

What is immediately seen is an xi ^2 dependence of the results, which change up to 100% 
even in this narrow kinematical range. Such a dependence is found at both the experimental 
and the model scale. The slope of the surface in the right panel of Fig. 2 is inverted w.r.t. that 
in the left panel. It is not a surprising feature, due to the different evolution properties of the 
numerator and denominator in Eq. (I12p and consistent with the evolution calculated in Ref. 
[To] . This X dependence, found at the hadronic scale as well as at high Q‘^, can be attributed to 
different dynamical and kinematical properties: 

1. both the numerator and the denominator vanish quickly with Xi through the valence region, 
but the latter vanishes faster, mainly due to the kinematic constraint xi + X 2 < 1 of the 
dPDF, a quantity appearing only in the denominator. The LF model correctly reproduce 
such a kinematical constrain; 

2. the correlations introduced by the LF dynamics and effective both in the Xi and k± de¬ 
pendence; 
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3. the correlations induced by the pQCD evolution in the valence region. 

In order to be more intuitive let us restrict to two different extreme scenarios: 

i) At very low-x gluons are strongly dominating (this is the hypothesis in [T7], partially 
corrected in m), so that it is enough to consider = g. Assuming, in addition, a 

fully factorized approach: Dgg{x, x' ,'k±) = Fg{x)Fg{x')g{]<-±), becomes: 


(Teff{xi,x\,X2,x'2) ^ (T^ff = \ ’ ( 1 ^) 

J r(k±)^ 

A similar assumption is used in Ref. m to obtain an estimate of a^fj which turns out 
to be about twice the experimental value. Obviously, the validity of Eq. lfH]) is spoiled by 
correlation effects and restricted to very low-x*. The problems related to the uncorrelated 
ansatz are discussed in a number of papers (see, e.g.. Ref. O [29l [30l [31] ). In particular, in 
the valence region this assumption is not supported by model calculations laiiKio] and it 
is certainly untrue in pQCD, being also spoiled by QCD evolution. In other words, several 
arguments lead to the conclusion that, in general, a^fj should be Xi dependent, namely: 
breaking of the factorization ansatz; the QCD evolution; contribution of more than one 
parton type (not only gluons as at very low xQ to the DPS cross section. 

ii) Let us now consider a simple way to reduce the results of our calculation to a fully factorized 
approach to dPDFs, following the hypothesis often assumed (cfr. Eqs. ([6|), ([7])): 


Duu{xi,X2-,'k±,g.o) = u{xi,g.o)u{x2,Ho)fuuO^±) , (15) 

where a natural definition for the “effective form factor”, /uu(kj_), in our approach, is 

fuui^x) = \j dxidx2Duuixi,X2-,'k±, Ho) , (16) 

a quantity which turns out to be scale independent. Within this approximation, Eq. (|12p 
yields: 


^ e//(^1 ) ) ^2 ) 3^2 )/ Tq ) ^ ^ e // 


81 


~ 10.9 m6, 


(17) 
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a value which turns out to be independent on the momentum scale Q and on the longitudi¬ 
nal momentum fractions Xi,x[. It is remarkable that it compares reasonably well with the 
sets of data shown in Fig.l. Of course the validity of our simplified result is restricted to 
the valence region where the model is predictive and the numerical estimate is connected 
to the ability of the model to capture (in its wave function) the correct average distribution 
of the valence quarks in transverse space. 

The x-dependence we are discussing does not emerge from the present data, probably not 
accurate enough. Our study points out therefore to an experimental scenario where more precise 
measurements in narrow Xi regions could shed new light on the structure of the proton and on 
the nature of hard proton-proton collisions. If the x-dependence is seen, one will gain, through 
(Teffi a first indication of double parton correlations and a fresh look at the 3D proton structure. 

4 Conclusions 

We have calculated the effective cross section (T^ff within a relativistic Poincare covariant quark 
model. Extracted from proton-proton scattering data by several experimental collaborations in 
the last 30 years, Ue// represents a tool to understand double parton scattering in a p-p colli¬ 
sion. Our investigation predicts a behavior of a^ff which, when averaged over the longitudinal 
momentum fractions Xj, is consistent with the present experimental scenario, in particular with 
the sets of data which include the valence region. However, at the same time, an x* dependence 
of <Xe// is found, a feature not easily read in the available data. We conclude that the measure¬ 
ment of CTeff in restricted Xj ranges woud lead to a first indication of double parton correlations 
in the proton, addressing a novel and interesting aspect of the 3D structure of the nucleon. 
The analysis of peculiar processes where these effects could be most easily seen, as well as the 
extension of the model to obtain a better description of the low-x region, presently studied at 
LHC, are in progress. 
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